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Abstract 

The coherent quantum dynamics of an electron in the quantum-dot ring 
structure under the resonant electromagnetic pulse is studied theoretically. A 
possibility of the selective electron transfer between any two dots is demon- 
strated. The transfer probability as a function of the pulse and dot parameters 
is calculated. It is shown that this probability can be close to unity. The fac- 
tors lowering the transfer probability in real systems are discussed. The results 
obtained may be used in the engineering of novel nanoelectronic devices for 
quantum bits processing. 

1 Introduction 

Rapid development of nanotechnology and unavoidable progress in miniaturization 
of the basic elements of contemporaneous microelectronics give rise to the new field 
of investigations, the physics of low- dimensional structures and nanoelectronics. In 
recent years, the great advances have been maid in the fabrication of nanostructures 
and in the study of their properties [1],[2]. Much attention is now paid to the 
quantum dots (QD) - "artificial atoms" [3] - which combine the properties of real 
atoms with the properties imposed by the process of fabrication. A possibility to 
use the QDs for quantum information processing is now extensively discussed [4]- 
[7]. In these hypothetical devices, known as quantum computers, the quantum 
information is encoded in ground and/or excited orbital states of electron in the 
QD or in electron spin degrees of freedom. In principle, it will be possible to realize 
the quantum algorithms in such systems. To achieve this goal, however, one has 
to overcome some challenges concerning with initialization, processing, readout and 
storage of quantum information. One of the main problems here is a coherent control 
of evolution of electron states under the influence of external fields. 

The behavior of low- dimensional objects is governed by quantum effects. Coher- 
ent evolution of one-electron states in a double-dot system upon the influence of a 
resonant laser pulse presents a good example of such phenomena. As was recently 
shown [8], [9], the pulse parameters (the frequency, duration, and amplitude) can be 
chosen so as to drive the electron, localized initially in the ground state in one of 
the dots, to the ground state of another dot via the delocalized state, common for 



both dots and used as the "transport state". If the states locahzed in different dots 
arc viewed as the Boolean states and 1, the electron transfer between them may 
be considered as the unitary operation NOT [8]- [10]. 

Recently, an attempt has been made [11,12] to generalize the results obtained 
for a double-dot system [8] to a chain-like multi-dot system. It was shown that it 
is rather difficult to implement the selective (addressed) electron transfer between 
arbitrary two QDs. This is because the probability amplitude to find an electron in 
a given QD depends strongly on the location of this QD in a linear chain of QDs 
with free boundary conditions. 

In this work, we demonstrate the possibility of selective electron transfer between 
arbitrary two QDs in the QD's ring structure under a resonant laser pulse, the local 
bias voltages being applied to those two dots. We derive an analytic expression for 
the transfer probability which takes into account the possible deviations of QD's 
and/or bias parameters from ideal ones as well as the detuning of the laser pulse 
from the resonance. 

2 The model 

We consider a system composed of N identical QDs arranged in a ring, see Fig.l. 
We suppose that there are at least two size-quantized levels in each QD. One of 
them, with the energy Si, corresponds to the ground state |1)„ localized in a given 
QD with the number n — 1; N. If the value of Si is close to the minimum of 
the potential energy of an electron in the QD, and the height and/or width of the 
potential barriers between QDs are large, then the ground state wave functions 
of neighboring QDs overlap weakly because of the strong localization of the wave 
functions (r|l)n in the corresponding QDs. In this case, the ground state of the 
whole system may be considered as A^-fold degenerate with respect to the electron 
localization in the QD system. We suppose that the excited levels |2)„ with the 
energy £2 (not necessarily second in the level numeration) in each of the QDs lie 
close to the barrier edge. In this case, the wave functions (r|2)„ of the neighboring 
QDs overlap strongly, resulting in the electron tunneling between QDs and splitting 
of the excited levels into the subband of N levels, each being delocalized over the 
QD system. Since in the following we will use the resonant (with respect to the 
external time-dependent field) approximation, we neglect all levels whose energies 
are far from ei and 62- 

The Hamiltonian of an excess electron added to the conduction band (over the 
filled valence band) of such a QD structure, has the form 

N N N 

n=l n=l n=l 

where a5^„(ai,„) and a2^„(a2,n) are, respectively, the operators of creation (annihi- 
lation) of an electron in the ground and excited states of an isolated QD with the 
number n; V is the matrix element of electron tunneling between the excited states 
of neighboring QDs. Note that a2,iv+i = 02,1 since the QD structure has the ring 
form. We don't show the spin indexes explicitly in the Hamiltonian (1) since we 
consider a one-electron problem. 

Let an electron be initially localized in the ground state |l)n,^ of a QD with the 
number rii. In the absence of an external field, the electron lifetime in the state 



\l)ni is very long because of weak overlap of the ground state wave functions of 
neighboring QDs. We assume this time to be longer than all other characteristic 
times of the problem. We wish to realize the selective electron transfer to the ground 
state 1 1)^2 of a QD with the number n2, i. e., to change the location of an electron 
in the QD system. 

To realize the selective electron transfer between the QDs with the numbers rii 

and n2, we apply equal local bias voltages to those QDs, thus shifting the energies of 
their ground and excited electron states (si —>■ Si+U, £2 S2+U). The Hamiltonian 
becomes 

H^Ho + C/(a+„^ai,„, + a+„^ai,„J + U{al^^a2,m + at^n^n^), (2) 

where, for the sake of simplicity, we ignore the changes in the matrix elements 
of electron tunneling from the QDs with the numbers rii and 712 to neighboring 
QDs. This approximation is valid if the local biases are small compared to V. For 
definiteness, we will consider U < and f/ < (in general, U ^ U, although they 
are of the same order). To diagonalize the Hamiltonian (2), we go to operators 
at = En=i Ck,na2,n = 1; -■; N). Wc have 

N N 

H ^£iJ2 «M«i,n + U (a+„^ai,„i + a+^^ai,^^ + ^katak, (3) 

n=l n=l 

where the energies Ek of the delocahzed levels and the coefficients Ck,n are deter- 
mined by the set of equations 

EkCk,n = ^lCk,n — y{Ck,n-l + Cfe,n+l) + UCk,n{5nni + ^nn^)-: k = 1] ...;N, (4) 

where Ck,N+i — Ck,i- The coefficients Ck^n are the probability amplitudes to find 
an electron occupying the k-th state in the excited state |2)„ centered in a QD with 
the number n. They satisfy the normalization condition J2n=i \Ck,n\'^ = 1 for any 
k^l;...;N. 

Expanding Ck,n into a Fourier series, Ck,n — ;^ C'fe,nexp(i27r mn/N), we 
obtain from Eq. (4) the following relationships between Ck,ni and Ck^n^' 

Ck,ni — ■^kCk,ni + BkCk,n2-> Ck,n2 — B^Ck^ni + ^kCk,n2^ {^) 

where 

^ ^ _^ 1 B ^--Y^ exp(i27r(ni - n2)m/N) 

^~ N l^^Ek -62 + 2Vcos{2nm/ NY ^~ N Ek - €2 + 2V cos{2nm/Ny 

(6) 

It follows from Eqs. (5) and (6) that |Cfc^„J = IC^^njI for any delocalized level k 
from the excited subband. This is important for the following consideration. Note 
that the coefficients Cfc,ni and Cfc,n2 rnay be chosen real, so that Cfc,ni = ±C'fc,n2- 

Let the ac electric field E(i) = EoCOs(Qi) be imposed on the QD system. The 
field frequency Q is close to the difference between the energy of one of the levels 
from the excited subband, Etr (in the following called as the "transport level"), and 
the ground state energies of an electron in the QDs with the numbers ni and n2 
(hereafter we set the Planck constant h — 1). In the resonant approximation [8,9], 



the Hamiltonian reads 



n=l n=l 
N , ^ \ f 1 ^ 1 

E (si + U{5nni + (^n^Jj ai;„ai,„ + Etratrtttr - \ 2 ^xp {-i9.t) E A„aii;ai,„ + /i.e. \ , 

n=l ^ ^ n=l J 

(7) 

where p„ = (ir|p|l)„ are the matrix elements of the momentum operator; A{t) is 
the vector potential (we use the Lorenz gauge with zero scalar potential and neglect 
the interaction term quadratic in the vector potential), m* is the electron effective 
mass. In Eq. (7), we introduced A„ = — ^^^EoPn making use of the well-known 
relationship between the vector potential and the strength of an electric field with 
the frequency ft and the amplitude Eq. 

Here we point on a relationship between An and the coefficients Ctr,n in the ex- 
pansion of the delocalized transport state \tr) = E^=i C'tr,n|2)n in the states |2)„. 
From the definitions of A„ and p„ one obtains 

ie ^ 

>^n = --T7,^oEC:r,n'-n'{2\mn. (8) 
n' 

Since the wave functions (r|2)n of the excited states of QDs are centered in the 
vicinity of the corresponding QDs, and the ground state wave functions are localized 
in the QDs, one may suppose that „/ (2| p = „ (2| p , then it follows from 

Eq. (8) that A„ = \Ctr,n, where A = -^^EqP and p = „ (2| p . 

We note that p 7^ (i. e., A 7^ 0) only if a certain relationship between the 
symmetries of the wave functions (r|l)n and (r|2)„ takes place. For example, p = 
if both those functions have s-symmetry, while p 7^ if one of them has s-symmetry 
and another has p-symmetry. Besides, for the value of A„ to be independent of n 
(this is needed to meet the equality |A„J = |A„2| which follows from the condition 
IC'ir.ml = |Ctr,n2l obtained earlier and to optimize the electron transfer between 
QDs), the vector p (not only its absolute value) should be independent of n. This 
is so if, e. g., the functions (r|l)n have s-symmetry and the functions (r|2)n have 
p^-symmetry, where z axis is perpendicular to the QD ring plane, see Fig.l. In this 
case, Eq should have non-zero component along z axis in order A 7^ 0. 

In the resonant approximation, the evolution of the electron wave function 

TV 

|*(t)) = ^ B„(t)|l)„exp[-i(£i + U{5nn, + Snn^M + Btr{t)\tr) exp{-lEtrt) (9) 
n=l 

is governed by the non-stationary Schrodinger equation 

Ant)) 



dt 



-H{t)mt)) (10) 



with the Hamiltonian (7) and the initial conditions -B„(0) = Snm, Btr{0) = 0. Our 
goal is to calculate Bn(t), Btr(t) and thus to find |^(t)). The probability to find an 
electron in the ground state of the QD with the number n at a time t is equal to 

Pn{t) = \Br,{t)\\ 

We choose the lowest level of the excited subband as a transport level. This choice 
is motivated by the following considerations. First, since the excited states of the 
QDs are close to the barrier edge, some of the upper levels of the excited subband 



may belong to the continuum spectrum. In contrast, the energy of the lowest excited 
level becomes smaller than 62, and hence the corresponding wave function remains 
localized in the QD system (although delocalized between different QDs). Second, 
the lowest excited level a.t U = is non-degenerate for arbitrary and remains 
non-degenerate at C/ 7^ 0, while other excited levels at U — constitute the set of 
the doubly degenerate states (except for the highest level at even N). AtU 7^ 0, this 
degeneracy is lifted, see Figs. 2 and 3. but the energy separations within the doublets 
are small compared with the doublet separations themselves, so that choosing any 
but the lowest of the excited delocalized levels as a transport level makes it difficult 
to tune the laser to the resonance. 



3 Results and discussion 

We define the resonant frequency and the detuning from the resonance as Q,. = 
Etr — (^1 + U) and 6 = Q ~ Qr, respectively. From Eqs. (7), (9), and (10) we obtain 
the set of the differential equations for the coefficients Bn{t) and Btrif): 



B„{t) = l^X*^Btrit) exp [l(^S-U{l- Snm - Snn2)) t 
Btr{t) T, KBn{t) exp y-i [5 -U{1- 5nm - ^nna)) t 



(11) 



where we took into account that the states \l)n and \tr) are the eigenstates of the 
stationary Schrodinger equation with the eigenvalues Si -\- U{5nni + ^nni) and Etr, 
respectively. 

Since the shift of the QD levels U caused by the local bias voltages applied to 
the QDs is finite and detuning from the resonance is small (ideally, 5 — 0), one 

can take \S\ « \U\. Moreover, we assume that the inequalities \S\ << |A| and 
|A| << \U\ are satisfied, so that \6\ << |A| << \U\. Then, as follows from Eq. (11), 
the characteristic time ~ 1/|A| it takes for the coefficients i?„^(t), Bn^it), and Bfr to 
vary is much longer than the corresponding time ~ ^/\U\ for the coefficients Bn{t) 

with n ^ ni,n2, and hence we have |-B„^ni,n2l ~ (I''^l/I^l)l-Sn=ni,n2l << \Bn=nun2\- 
Thus, it is sufficient to restrict ourselves with n — rii and n = ^2 in the sum (11): 

Bn, (i) = i^Bt, (t) ■ exp it5t) , 

B^At)=^^Btrit)-exp{i6t), (12) 
Btr (t) = i^Bn^ (t) ■ exp {-iSt) + i^B^^ (t) ■ exp {-iSt) , 

Bn^n„n,{t) = l^Btrit) C^K^ " U)t). (13) 

In this manner, the (A^+l)-level problem can be reduced to the 3-level problem since 
there are only three levels (|l)ni, |l)n2) and \tr)) relevant for the quantum dynamics 
of the system at the resonance or close to it. This problem was recently solved by 
us for the case 5 7^ 0, |AnJ = jAnjI in Ref. [8] and for the case 5^0, jA^J 7^ |A„2| 
in Ref. [9]. The results obtained in Refs. [8,9] for the probabilities p„(t) can be 
directly applied to the problem treated here. One can see [9] that at 5 = 0, the 
probability of electron transfer between QDs is 

PnM = fn^%#4Vsin^Kt), (14) 



where cur — y^jA^|j^~-|-]A^^/4, and hence the selective electron transfer between 
QDs takes place in a time T = t:/2ujr if |AnJ = |A„2|- A deviation of 5 from zero 
and |A„J from |A„2| causes the value of Pn2{t) to decrease. 

In this work, we account for a possible differences in the voltage biases applied 
to two selected QDs. We take Un^ ^ and 17^ 7^ Un2: then Eqs. (12) and(13) 
become 



Bn2 (t) 
Btr {t) : 



i^Btr (t) ■ exp {i6t) , 
i^Btr (t) ■ exp {i{6 - Ae)t) , 

i^B^^ (t) ■ exp {-i6t) + i^Bn^ (t) ■ exp {-i{6 - Ae)t) , 
A* 

Bn^nuniit) = i^Btr{t) exp(i((5 - t/„Jt), 



where now VLr = -^tr — (^i + f^jj and we designated Ae = C7, 



n2 • 



Using the following substitutions, 



B. 



ni 



Br^2 (t) 



(t) = Br,, (t) ■ exp (ift) , 
(t) ■ exp [-i^t 



B. 



Btr (t) = Btr (t) ■ exp [-i (S -f)t), 



we have from Eq. (15): 



(15) 



(16) 



(17) 



Bn, (t) 



• Ae 



Kit) 



l^Btr (t) 



it)-tfBn,it)=l^Btr (t) 

[ Btr{t)-i{5-Ae/2)Btr{t) 



{t)+l-- 



2 ^n2 



it). 



(18) 



Next, we express Bn,{t) in terms of Bn^{t), Btr{t), and Btr (t) and Btr{t) in terms of 

Bn^{t) and [t) from the third and second equations of the set (18), respectively. 
Inserting these expressions into the first equation of the set (18), we obtain the 
equation for the coefficient Bn^{t): 



Bn,{t)-z{6-Ae/2)Bn,{t) + 



Ae 
2 



[Anil -[Anal ((5-A£/2)A£ 
4 2 



s„2(0 = o. 



(19) 



Taking into account that -B„^(0) = 1 and -B„2(0) = -Btr(O) = 0, we have the following 

initial conditions: Br,^(0) = -B„2(0) = 0, B^^iO) = -A„iA*2/4. 

The solution of Eq. (19) can be found exactly. It is, however, too combersomc. so 
here we restrict ourselves to the approximate solution for the case that the deviations 
from the ideal case are small, i. e., |5|/|A|, |A£|/|A|, |AA|/|A| << 1, where AA = 

\^ni I |An2 |- 



Br,, (t) = r^^^ exp -t ,171 

• {1 - [cos {2ujRt) - sin {2ojRt)] exp (z|)} 



(20) 



|2 



where 5 = 5-^ + fA^p^Uf^ and ur = y\\n,f + \\n/ + P + {A6f. Then 



we have for the probabihty of the electron transfer: 



Pn, (t) = \Bn, 



2 



2 


1 "^^1 




Ana 




1 '^^1 1 


2 

+ 


A 


12 


+ 


;Ae)2 



sin4(a;Ri) + sin^ (f ) cos{2ujRt) + ^ sm\2ujRt) - ^ sin (f ) sm{2ujRt) 

' (21) 

For |5|/|A|, |A£|/|A|, |AA|/|A| << 1, the maximum value of is 

^ 1 - [(AA)^ + (As)^ + (ttVs) (5 - A£/2)^]/|A|^ (22) 

at T = 'n/2u}R, in agreement with Ref. [9]. 

If at f = an electron is placed in the superposition of states and |l)n2) i- 6., 
= a and B^M = (so that |q;|' + = 1, i. e., Bn^ri,,nM = ^tr(O) = 0), 
then for the ideal structure (Ae = AA = 0) and at exact resonance {5 = 0) one 
obtains Bn^{T) = —f3 and Bn^lT) = —a. Thus, if the electron states and 
\l)n2 localized in different QDs are considered as the Boolean states |0) and |1), 
respectively, so that their linear combination is viewed as a quantum bit (qubit), 
then the resonance pulse with duration T is equivalent to the unitary quantum 
operation NOT combined with the simultaneous change of the qubit phase by tt: 



UidealiT) 1^(0)) = U^deal{T) K = |*(T)) = K exp(z7r). (23) 





At non-zero values of 6, Ae, and AA, the fidelity of this operation, F = (UidealiT)'^ (0) U (T)^(O) 

\-p*Br,,{T) - a*Bn^{T)f, differs from unity. 

For the state-of-the art nanotechnology, it seems very difficult to fabricate an 
ordered nanostructure composed of a great number of nearly identical QDs. It is 
worth to mention another physical system, for which the realization of the proposed 
scheme may be possible, an array of phosphorous donor atoms embedded in a sil- 
icon host [13], [14]. The modern techniques of the controlled implantation of the 
phosphorous atoms into a silicon substrate [14] allows, in principle, to fabricate the 
structures where the donors serve as the centers of electron localization. Unlike an 
"artificial" (QD) atoms, all "natural" atoms are identical, while possible differences 
caused by implantation defects can be minimized, by appropriate annealing. If all 
but one of donors are ionized, the one-electron model studied in this paper may 
be used to describe the evolution of an electron state. 

Finally, since an electron resides in a solid rather than in a free space, its unavoid- 
able interactions with other degrees of freedom can destroy the unitary electron evo- 
lution under external pulse. In particular, the processes of electron relaxation and 
dephasing result to the decoherence. This imposes some technological restrictions 
on the structure and material parameters [15]. A detailed discussion of decoherence 
effects in such structures will be presented elsewhere. 



4 Conclusions 



In this work, we have studied the unitary evolution of one-clcctron states in an array 
of coupled QDs. We have shown that the selective electron transfer between two 



arbitrary QDs under the influence of tlie resonant pulse is possible. The probabilities 
to find an electron in the ground states of those QDs as functions of time and the 
structure and pulse parameters were calculated. We have shown that the probability 
of electron transfer to the ground state of another QD may be very close to unity. If 
we consider the states localized in two QDs as logical variables "0" and "1", then the 
resonant pulse applied to an arbitrary linear superposition of those states (qubit) 
implements the quantum NOT operation with the simultaneous change of the qubit 
phase by vr. 

We believe that since the selective electron transfer by virtue of the "transport" 
state delocalized over all QDs may be performed as a "one-step" operation, it has 
an advantage over the scheme based on sequential " step- by-step" electron jumps 
between neighboring QDs, as was recently suggested in Ref. [16]. The results 
obtained in our work can be directly extended to other nano structures composed 
of a large number of centers of electron localization, e. g., the phosphorous donor 
in a silicon host, the atoms absorbed on a surface, etc. 

We are grateful to K.A. Valiev for his interest in this study and to S.A. Dubovis 
who participated this work at its early stage. 
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Fig. 1. The quasi-one-dimensional nanostructure composed of identical QDs 
arranged in a ring. The equal bias voltages are applied to the QDs with the numbers 
rii and n2. The selective electron transfer between those QDs is performed under a 
resonant laser pulse. 
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Fig. 2. The energy diagram of the A'^-dot structure shown in Fig. 1 shows the 
energies Ej. of the delocahzed states in the excited subband versus the local bias 
voltages U applied to the QDs with the numbers rii and n2- = 9, ni = 1, ^2 = 4. 
The solid and dash lines correspond to the non-degenerate and degenerate at t/ = 
states, respectively. 
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Fig. 3. The same as in Fig. 2 for N — 10, ni = 1, ^2 = 5. 



